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ABSTRACT  (UNCLASSIFIED) 


In  azimuth  processing  of  SAR  data,  the  slant  range  histories  of  the  targets  are  requited.  After  range 
compression  the  maxima  of  the  responses  of  a  point  target  are  on  the  so  called  range  migration  locus.  For 
spacebome  SAR  the  locus  traverses  a  number  of  range  bins.  This  locus  determines  where  to  get  the  samples 
needed  for  the  azimuth  compression  and  the  phase  function  which  must  be  used.  For  satellite  SAR  the  locus 
can  be  described  by  two  parameters,  the  Doppler  centroid  and  the  frequency  rate  These  parameters 
can  be  estimated  by  using  accurate  spacecraft  ancillary  data  and  by  utilizing  the  received  radar  data  alone. 
This  report  describes  the  effects  of  errors  in  the  parameters  used  by  the  processor  and  discusses  algorithms  to 
estimate  the  parameters.  Exan^les  of  results  from  processing  Seasat  data  are  shown. 

The  report  is  one  of  a  series  for  a  SAR  processor  for  the  ERS-1  data,  phase  2,  precision  processor.  The  work 
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SAMENVATTING  (ONGERUBRICEERD) 

Voor  de  azimut  compressie  van  SAR  data  moet  bet  verloop  van  de  afstanden  van  de  doelen  tot  de  satelliet  in 
de  tijd  bekend  zijn.  Na  bet  comprimeren  in  range  liggen  de  nuuuma  van  een  puntdoel  responsie  op  de 
zogenaamde  range  migratie  locus.  Voor  satelliet  SAR  doorloopt  deze  locus  een  aantal  afstandscellen.  Deze 
locus  bepaalt  waar  de  bemonsteiingen  voor  de  azimut  compressie  liggen  en  de  fase  functie  die  gebruikt  moet 
worden.  Voor  satelliet  SAR  zijn  twee  parameters  voldoende  om  de  locus  te  bescbrijven.  Deze  parameters 
zijn,  de  Doppler  centroide  fpc  veranderingssnelheid  van  de  Doppler  frequence  Deze  parameters 
kunnen  worden  geschat  uit  nauwkeuiige  hulpdata  van  de  satelliet  en  uit  de  ontvangen  radatdata  zelf.  Dit 
rapport  beschrijft  het  effect  van  fouten  in  de  parameters  die  door  de  processor  gebruikt  worden  en  bespieekt 
algoritmen  voor  het  scbatten  van  de  parameters.  Voorbeelden  van  resultaten  'veikregen  met  het  processen 
van  Seasat  data  worden  getoond. 

Het  rapport  is  £dn  uit  een  serie  voor  een  precisie  SAR  processm-  voor  ERS-1  data.  Het  bescbieven  weik  is 
gesubsidieerd  door  NIVR  (NIVR,  NRT  nr.  2802FE,  SAR  processor). 
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1  INTRODUCTION 

In  the  azimuth  processing  of  the  SAR  data,  the  slant  range  histories  of  the  targets  are  required. 
The  changing  slant  range  between  the  radar  and  the  target  as  the  radar  moves  past  is  r(r).  First,  the 
azimuth  samples  used  in  azimuth  compression  have  to  be  taken  from  different  range  gates  in  each 
pulse  as  dictated  by  r(t).  The  slant  range  also  determines  the  phase  history 


<>(r)  =  (4jiA)r(f)  (1.1) 

where  X  is  the  radar  wavelength. 

The  phase  is  used  in  generating  the  azimuth  matched  filter.  Accurate  determination  of  the 
functions  is  required.  Errors  in  these  functions  degrade  the  image  quality.  The  functions  can  be 
estimated  by  using  accurate  spacecraft  ancillary  data  and  by  utilizing  the  radar  echo  data  alone. 
Depending  on  the  scene,  with  utilization  of  the  raw  radar  data,  greater  accuracies  can  be  achieved. 
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2  THE  PHASE  FUNCTION  FOR  A  SPACEBORNESAR 

For  spacebome  SAR's,  the  random  motion  of  the  spacecraft,  which  induce  variations  in  die  {diase 
histories  of  the  targets  from  those  expected  in  a  smooth  flight  path,  is  usually  negligible.  So, 
motion  compensation  (the  techniques  for  removing  the  effects  of  the  phase  variations),  is  not 
required.  It  can  be  shown  [Barber  85],[lj  8S]  that  the  expected  phase  history  can  be  approximated 
by  a  function  with  two  parameters  c,  and  Cj 

(KO'Co+Cit  +  Cjt^  (1.2) 

where  t-Q  corresponds  with  the  moment  where  the  target  is  in  the  antenna  boresigbt  direction. 
The  instantaneous  Doppler  firequency  is  given  by 

/(0  =  -  ^  (1.3) 

2Tt  dt 

and 

fit)  ^  foe  fpKi  (1.4) 

The  azimuth  Doppler  history  is  thus  a  linear  chirp  centred  at  having  a  Doppler  rate  The 
two  terms,  fpc  and  fpg,  will  be  referred  to  as  the  Doppler  centroid  and  the  Doppler  fiequency 
rate,  reflectively. 

Thus,  the  range-migration  curve  r(t)  aad  the  jdiase  history  curve  ^r)  can  be  estimated  by 
determining  the  two  Doppler  parameters. 

In  the  sequel  we  use  the  following  ootations  and  simple  relations 
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3  THE  EFFECT  OF  PHASE  ERRORS 


The  accuracy  requirements  for  the  Doppler  parameters  are  determined  by  the  effects  on  image 
quality  resulting  from  using  erroneous  parameters  in  the  processor. 

3. 1  Doppler  centroid 

The  response  of  a  point  target,  in  the  range  compressed  domein,  is  a  linear  FM-signal 


s,(/)  =  j(/)  exp(/(K0)  (1.10) 

where  a{t)  is  caused  by  the  antenna  pattern.  The  antenna  pattern  is  by  construction  a  symmetric 


-PRF/2  < - Pwcetied  Bamlwidt  (PBW)  ■  ">  FltF/2 

beqaacjr 

Figure  1.1:  Frequency  ipectra  of  poini  target  response  and  azimuth  lefeience  ftmctioa 

one.  Assuming  a  target  with  no  variationf  of  bsKkacatter  as  a  function  of  aTmuitb  angle,  we  find 
that  a(0  is  a  symmetric  ftmctiaa.  Thus  the  ^ectrum  of  Ae  linear  FM-pnlse  weighted  by  Ae 
symmetric  function  a(0  is  syimnmric. 
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-PRF/2  < - Processed  Bandwidth  (PBW) - >  PRF/2 

frequency 

Figure  1.2:  Spectra  of  signal,  ambiguities  and  azimuth  reference  function 


The  azimuth  replica  is  a  linear  FM-pulse  with  constant  amplitude  so  the  azimuth  transfer  ftmction 
is  symmetric  and  flat  having  a  fixed  Processed  Bandwidth  (PBW).  We  get  a  maximum  signal 
output  when  the  reference  Doppler  centroid  is  equal  to  the  signal  centroid  /^.  By  using  an 
incorrect  Doppler  centroid  we  lose  signal  energy,  as  is  evident  from  figure  1.1  (see  also  A12  in 
Appendix  A).  The  white  noise  contribution  to  the  conelator  output  does  not  depend  on  the 
centroid  value  fpp  used.  Thus,  Dof^ler  centroid  errors  lead  to  a  decrease  in  the  signal-to-noise 
ratio  (SNR). 

The  actual  received  radar  sigrud  is  a  stochastic  signal  r(r,a),  consisting  of  the  superposition  of 
point  target  responses,  each  having  its  own  position  ('',.,0^,  amplitude  and  phase.  Note  that  we 
have  replaced  time  r  by  azimuth  a,  which  is  proportional  to  t. 

Now,  consider  the  one-dimensional  signal  s^a)ms{r4),  where  rsconstuit.  We  get  an  estimate  of 
the  PSD  of  this  signal  by  taking  its  Fourier  transfmm  and  determining  S^CflS  Xf)-  A  better 
estinuue  is  obtained  by  averaging  a  number  of  spectra  adjacent  in  range. 

Figure  1.3  shows  three  average  image  ^reclra  (qrectra  of  azimuth  cott^nessed  Seasu  data), 
corresponding  with  three  different  Doppler  centroid  values  used  by  die  processor.  The  qiecttum 
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in  the  middle  corresponds  with  the  situation  where  die  Doppler  centroid  value  used  by  the 
processor  is  the  Doppler  centroid  value  of  the  signal  (/£>7>=/^)-  We  see  that  the  image  spectrum 
is  symmetric.  The  upper-  and  lower  spectra  correspond  with  Doppler  centroid  values  fpp  * 

The  matched  filter  used  by  the  processor  has  a  fiequency  spectrum  that  does  not  centre  on  the 
maximum  return  of  the  antenna's  azimuth  power  pattern  which  results  in  a  asymmetric  image 
spectrum. 

The  PRP  is  in  the  order  of  the  frequency  range  traced  out  by  the  Doppler  of  the  target  during  its 
traversal  through  the  antenna  main  lobe.  So  there  will  be  aliased  energies  in  the  processed 
bandwidth  (PWB)  from  the  ambiguities.  The  minimum  amount  of  ambiguity  energy  is  in  the 
PBW  when  the  processor  Doppler  centroid  fpp  equals  the  correct  Doppler  centroid  fpp.  Thus 
Doppler  centroid  errors  lead  to  a  decrease  in  the  signal-to-azimuth  ambiguity  level  ratio,  as  is 
evident  from  figure  1.2.  The  azimuth  ambiguity  and  a  definition  of  a  signal-to-ambiguity  ratio 
(SALR)  is  discussed  in  [Tomiyasu  80], 

In  [Li  85]  (see  also  [Curlander  82]  the  nominal  orbital  and  radar  system  parameters  for  the  Seasat 
S  AR  have  been  used  to  generate  curves  fat  fpp  and  fpp  at  the  radar  boresight  as  a  function  of  the 
true  anomaly  of  the  orbit  for  the  cases  with  no  spacecraft  attitude  errors  and  with  1*  eiror  in  yaw 
and  r  error  in  pitch.  These  attitude  errors  can  cause  Doppler  centroid  shifts  of  more  than  lOOOHz 
(the  PRF  is  -  1647Hz).  Also  evaluated  and  plotted  were  the  losses  in  SNR  and  SALR  as  a 
function  of  the  Doppler  centroid  error.  The  plots  show  that  the  SALR  deteriorates  much  faster,  as 
a  function  of  Doppler  centroid  error,  than  the  SNR.  For  losses  in  SALR  and  SNR  both  less  than 
3db,  the  Doppler  centroid  error  must  be  less  than  -150Hz,  which  corresponds  with  a  spacecraft 
attitude  uncertainty  of  less  than  -0.2".  For  the  Seasat  experiment,  the  attitude  is  generally 
determined  to  better  than  -0.1*.  There  are  occasions,  however,  when,  the  attitude  is  not  known  to 
better  than  -0.7"  due  to  malfunctioning  of  the  attitude  determination  system,  the  attitude  is  not 
known  to  better  than  —0.7".  In  diese  cases  use  of  the  resulting  erroneous  estimates  causes 
significant  losses  in  SNR  and  SALR.  Furthermore,  there  are  applications  using  SAR  data  that 
require  very  accurate  estimates  of  fpp.  These  accurate  estimates  can  be  obtained  using  the  radar 
echo  data  themselves. 

In  Appendix  A,  we  have  shown  that  the  point  target  response  is  shifted  in  time  (azimuth)  and 
range  fiom  its  correct  position  when  using  an  incorrect  Doppler  centroid  value. 
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The  shift  in  the  azimuth  dimension  is 


^  _  *1,  -  aj,  for  -  foe 
2^2/  foK 

and  the  shift  in  the  range  dimension  is 


(1.11) 


a|,-aj,  ^  -X 
~2 


(Jar  +  foe) 


(fof  "  foe) 


(1.12) 


The  shift  in  the  azimuth  dimension  is  illustrated  by  Figure  1.4.  Shown  are  8  images  of  the  Seasat 
scene  containing  the  Goldstone  antenna.  For  each  image  a  different  processor  Doppler  centroid 
value  fup  was  used.  The  frequency  values  between  two  adjacent  images  differ  by  17.0  Hz.  For 
each  image,  range  goes  from  left  to  right  and  the  azimuth  is  in  the  vertical  direction. 

3 . 2  Doppler  frequency  rate  fpg 

Consider  a  point  target  which  at  r=0  is  in  the  antenna  boresigbt  direction.  To  image  the  target  we 
can  use  any  interval  of  the  target  locus  within  the  Processed  Bandwidth  (PBW),  e.g.  the  interval 
with  length  and  centre  x  (see  Figure  1.5).  The  azimuth  resolution  is  determined  by  the  aperture 
(TJ.  If  an  incorrect  Doppler  firequerrey  rate  is  used,  the  maximum  target  response  will  occur  at 
azimuth  (time)  and  range  not  corresponding  with  the  target's  azimuth  and  time.  This  is  the  case 
because  the  (wrong)  reference  locus  will  approximate  the  target  locus  best  at  a  part  different  from 
the  one  sketched  in  Figure  1.5.  The  reference  locus  does  not,  however,  completely  coincide  with 
the  target  locus.  So  (slightly)  wrong  data  is  taken  and  a  (slightly )wrong  phase  reference  frinctioa 
is  used.  The  resulting  point  target  image  is  shifted  from  the  true  position,  the  response  is 
broadened  and  the  peak  gain  is  reduced.  These  effects  are  functions  of  the  cumulative  phase 
errOTS.  The  shift  of  the  point  target  response  depends  on  t  (the  azimuth  centre  value  of  the 
processing  interval). 

The  effect  of  an  error  in  the  Doppler  frequency  rate  in  a  n-look  SAR  processor  is  that  each  target 
will  appear  as  n  separate  targets  in  the  image. 

The  shift  in  time  is  (see  [Li  85]) 
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-  fott-frr  5/o*  ,, 

5/= - T  =  - t  (M3) 

frr  frf 

where  is  the  correct  Doppler  frequency  rate,  f^p  the  frequency  rate  used  by  the  processor  , 
^fop  the  rate  error  and  x  the  look's  centre  time. 

For  4-looks,  we  find  a  time  shift  of  the  two  images  corresponding  with  the  extreme  looks  LL  and 
RR  (see  figure  1.5)  of 


frP  ^ 


(1.14) 


where  r„  is  the  maximum  available  azimuth  integration  time. 

Using  typical  Seasat  values,  we  get  for  the  displacement  between  two  look-images  Sx  =v6t 


-  25  5/o»  m 


(M5) 


» 
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range  r 


where  =  2,48  sec,  v  -  7000in/sec  and  f^p  =  ■=  530Hz/sec. 

The  resolution  of  the  four  look  Seasat  image  is  -25m.  If  we  constrain  the  image  displacement  to 
less  than  8m  than  ifop  <  0.3  Hz/sec.  The  corresponding  cumulative  phase  error  in  each 
subaperture  is  -5*.  This  phase  error  generates  acceptable  degradation  in  each  lotdt's  point  target 
response. 

The  Doppler  rate  fpg  is  not  as  sensitive  to  the  ^>acecraft  attitude  fluctuations  [Curlander  82]  as 
the  Doppler  centroid.  An  analysis  of  the  euor  sources  in  the  estiination  of  fpg  using  this 
technique  indicates  that  uncertainty  in  the  velocity  of  the  satellite,  is  the  major  fpg 
estimation  error  source.  It  has  been  determined  that  fcr  Seasat  a  satisfactoty  value  for  fpg  is 
within  0.75Hz/sec  of  the  true  fpg.  This  accuracy  level  can  be  obtained  if  v,  is  known  to  better 
than  0.5m/sec  in  each  conqranent.  Thus  using  the  ancillary  data  approach  may  provide  an 
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adequate  esdmats  for  fpg.  For  future  real  time  SAR  processors,  when  only  approximate  estimates 
of  the  ancillary  spacecraft  data  are  available,  the  expectation  is  that  such  accuracies  will  be 
difficult  to  obtain.  Use  of  an  autofocus  method,  therefore,  can  be  advantageous. 
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4  DOPPLER  PARAMETER  DETERMINATION  TECHNIQUES 

Report  [Otten  89]  presents  a  study  of  available  literature  on  phase  emv  reduction  techniques, 
supplemented  by  comments  of  general  nature.  For  satellite  SAR  the  i^ase  history  can  be 
approximated  by  a  function  with  the  two  Doppler  parameters:  the  Doppler  centroid  and  the 
Doppler  rate  fiy/g.  The  two  techniques,  to  estimate  these  parameters,  described  in  this  report 
achieve  accuracies  which  can  exceed  the  system  performance  requirements.  The  methods  to 
estimate  the  two  Doppler  parameters  and  are  referred  to  as  "clutterlock"  and  "autofocus" 
respectively.  The  methods  are  computationally  more  demanding  than  the  parameter  estimation 
based  on  accurate  spacecraft  ancillary  data.  So,  when  the  potentially  more  accurate  estimates  are 
not  needed,  we  use  the  ancillary  data.  Moieover,  both  the  clutterlock  and  autofocus  methods  need 
initial  estimates.  These  initial  estimates  are  ^vided  by  the  values  determined  from  the  ancillary 
data  or,  in  cases  where  such  data  are  not  available,  by  estimates  from  simpler  implementations  of 
the  clutterlock  and  autofocus  techniques. 

4. 1  Doppler  parameter  estimation  derived  from  spacecraft  ancillary  data 

Figure  1.6  shows  the  imaging  geometiy  involved  for  a  spacecraft  SAR.  The  range  from  satellite  to 
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SATELLITE 


Figure  1 .6:  Geomefiy  of  relative  position  between  sensor  and  point  taiget 


target  is  given  by 


HO=rj{t)-  rs(0  (1.16) 

whereat)  refers  to  the  position  vector  of  the  target  on  the  Earth's  surface  and  r^t)  refers  to  (he 
position  vector  for  the  satellite. 

The  target's  motion  can  be  modeled  as  a  point  on  the  surface  of  a  rigidly  rotating  Earth,  so 


Ho  “  ^0)  +  vj<0)(  +  i 


(1.17) 
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where  ^  (0)  =  m,  xiy  (0),  is  the  Earth's  angular  velocity  vector  and  Oj  (0)  is  the  acceloation 
vector  of  the  target. 

For  small  r,  the  satellite  motion  is  modeled  as 


r^t)  =  ^0)  +  V5{0)r  +  i 


(1.18) 


where  U5<0)  and  35(0)  are  the  velocity  and  acceleraboo  vectors  of  the  satellite. 
Substimting  1.17  and  1.18  into  1.16 


r(t)  =  F(0)  + 


?(0)f  -  i  a(Q)fi 
2 


(1.19) 


where 

K0)=  ^0)-  ^(0) 
u(0)=  u^O)-  05(0) 

3(0)=  3^0)-  3,(0) 

Inserting  1.19  into  r(f)  =  V  rtt)  ■  rtf)  and  retaining  terms  to  fl  we  get 


)f^ 


(0)  o(0)|f  + [0-^(0)  -  r(0)  a(0)]f^ 


f(f)  =  f(0).^ 

the  second  square  root  term  is  small  compared  to  1,  so  we  can  use 


f^O) 


We  find 
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r(f)  =  r(0)  + 


r(0)  5(0)  u^(0)-F(0)  a(0), 


r(0) 


-/  + 


2r(0) 


(1.20) 


Using  l.S  through  1.9  we  get 


2a,  ^  2  r (0)  5(0) 
A.  X  r(0) 


fox  ~ 


X 


2  \)^(0)-r(0).a(0) 

X  r(0) 


(1.22) 


The  spacecraft  position  ^0)  and  velocity  Vs(.0)  are  typically  known  from  the  orbit  tracking  data 
and  the  acceleration  is  easily  determined  assuming  negligible  drag  (i.e.  purely  gravitational 
acceleraticn). 

The  target  location  r|<0)  is  found  by  solving  the  vector  equation:  ^0)  =  ^(0)  .f.  riO)  Ea  assuming 
some  earth  model  where  is  the  boresight  unit  vector  and  t(0)  is  the  slant  range  from 
satellite  to  target.  The  position  of  the  target  on  the  eaith's  model  and  the  rotational  velocity  vector 
(0,  determine  the  velocity  vector  v,<0).  Since  both  the  satellite  and  target  postion,  velocity  and 
acceleration  are  known,  the  Doppler  parameters  and  therefore  the  range  migration  locus  and 
azimuth  refeience  function  can  be  determined.  Since  for  a  typical  SAR  the  antenna  is  rigidly 
mounted  to  the  satellite  platform,  the  attitude  variation  is  directly  translated  to  die  boresi^t  of  the 
radar.  Thus  must  be  adjusted  for  the  attitude  of  the  satellite.  The  accuracy  of  this  vectOT  is 
therefore  directly  d^ndent  on  the  uncertainty  in  the  qtacecraft  attitude  determination  system.  As 
mentioned  before  the  ancillary  data  tqiproacb  for  determination  of  is  unsatisfactory  uhereas 
it  may  provide  an  adequate  estimate  for  fpg. 

4.2  Doppler  centroid  deteimination,  CLUTTERLOCK 

The  technique  to  determine  dre  Doppler  centroid  bom  die  radar  echo  dida  is  baaed  on  the  fact  that 
the  radar  echo  dato  consisu  of  the  superposition  of  a  large  number  of  scatteren.  A  point  target 
response  in  die  azimudi  direction  is  a  linear  FM-sigiul,  centred  on  the  Doppkr  cemrmd, 
modulated  by  die  antenna  pattern.  This  antenna  pattern  is  by  design  symmetric  about  its  peak.  So, 
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if  in  a  scene  a  very  dcnninant  point  target  was  present  we  could  find  the  frequency  which 
corresponds  with  the  peak  response.In  reality  the  radar  signal  consists  of  a  siqierposititm  of  a 
large  number  of  scatterers.  The  resulting  signal  will  be  a  sample  function  from  a  randcnn  process. 
Now,  assume  a  target  area  which  has  a  homogeneous  scatter  density,  and  diat  the  scatterers  show 
no  variations  of  backscatter  cross  sections  as  a  function  of  azimuth  angle.  Then,  the  power 
spectral  density  (PSD)  of  the  azimuth  spectra  has  the  shape  of  the  antenna  power  pattern  and 
again  the  frequency  corresponding  with  the  peak  density  is  the  Dt^ler  centroid.  The  technique 
used  to  estimate  the  Doppler  centroid  consists  of  estimating  the  PSD  by  conqniting  the 
periodogram  of  an  azimuth  line  (that  is  at  constant  range).  To  obtain  statistically  consistent 
periodogram  spectra  a  number  of  periodograms  from  adjacent  ranges  are  averaged.  The  number 
must  not  be  too  large  since  is  a  function  of  range.  In  Figure  1.7  we  have  indicated  the 
domains  were  the  PSD  can  be  determined. 


n(a,r) 

©  ©  © 


Figure  1 .7:  SAR  encoding  and  processing;  1, 2  and  3  indicate  domains  from  which  the  data  can  be  used 
for  estimation  of  the  Doppler  centroid. 

Two  problems  have  to  be  treated  in  the  SAR  modelling.  One  is  to  relate  the  scene  reflectivity  and 
the  complex  ouqnil  image.  Secondly,  a  model  for  die  radar  scene  is  required.  These  problems 
have  been  treated,  extensive  references  can  be  found  in  [Madsen  89].  In  Figure  1.7,  Ojiajr)  is  die 
complex  radar  cross  section  lo^a/')F  w  ofa/’}  in  azimuth  and  slant  range  coordinates,  /,  ««d 
are  frequencies  in  die  azhnudi  and  range  dimeiisions,  H,  represents  the  encoding  of  the  scene, 
n{a/)  indicaies  thermal  noise,  /f^  is  die  range  compression  and  H,  is  die  azimi^  coo^iiession. 
For  the  PSD  estimation  we  use  windows.  There  will  be  a  lot  of  target  responses  dial  are  only 
partially  covered  by  die  window.  The  influence  of  this  fact  on  the  esdmatiaa  quality  is  more 
severe  as  a  target  response  efiecdvely  covers  a  large  area.  The  roost  confined  re^wnse  occurs  in 
the  final  domain,  die  image  domain  (point  3).  So  the  best  results  are  obtained  by  estimation  based 
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on  the  range-  and  azimuth  compressed  data,  which,  however,  requires  that  fpp  is  taken  into 
consideration,  leading  to  an  iteration  procedure.  The  main  advantage  of  using  die  raw  data,  point 
1,  is  that  the  estimation  jn-ocedure  will  not  depend  on  the  processor  setting  Estimation,  using 
data  from  this  domain  is  therefore  faster  than  by  utilization  of  image  data  and  can  be  used  in 
realtime  processing.  The  raw  data  can  also  serve  for  an  initial  estimate.  Since  we  are  concerned 
with  the  precision  processor,  we  have  chosen  for  the  estimation  of  the  Doppln  centroid  f^p, 
based  on  the  complex  image  data  (point  3). 

The  implemented  Dr^pler  centroid  estimation  method  is  the  so-called  A  E-method.  The  A  £- 
method  assumes  a  PSD  showing  symmetry  about  the  true  Doppler  centroid.  If  the  processor 
centroid  f^p,  equals  this  true  centroid  we  have  this  symmetry.  Then,  conrqiuting  the  energy  E/  in 
the  interval  at  the  lefthand  side  of  f^fp,  and  the  energy  E,  in  the  righthand  interval,  we  find  Ep^Ef. 
If  an  incorrect  processor  centroid  value  is  used,  the  energies  will  differ  A  £  =  A  E^^  ^  0. 

Using  a  number  of  different  processco'  centroid  values  we  find  a  curve  which  is  highly' line!ar.  By 
fitting  a  linear  function  of  frequency  to  the  conqiuted  A  £  values  we  can  determine  the  intercept 
with  the /£,y>  axis  (i.e.  the  frequency  value  where  A£  =  0).  This  frequency  value is  the 
estimate  of  the  Doppler  centroid 

Figure  1.8  shows  the  blockdiagram  of  the  Doppler  centroid  line  (the  Doppler  centroid  is  a 
function  of  range)  estimator.  The  clutterlock  algorithm,  called  the  A  £-method  operates  as 
follows: 

1 .  Based  on  a  quick  look  processing,  a  portion  of  the  raw  SAR  data  to  be  processed  into  an 
image  is  selected.  The  portion  is  divided  into  intervals,  the  data  in  an  interval  centred  on 
range  r  is  used  to  estimate  the  Doppler  centroid  /’'pcit).  Ova  the  relatively  small  area  used 
to  estimate  the  centroid  line,  the  dieoretical  variation  of  the  Doppler  centroid  as  a  function  of 
range  is  approximately  linear.  An  inioval  most  be  large  enough  to  be  able  to  perform 
statistical  averaging  to  get  stable  estimates  but  not  too  large  since  the  Doppler  centroid  may 
be  considered  constant  only  for  ^ort  range  intervals.  We  have  used  an  image  covering  an 
area  •>  16km  by  8km  (1024range  x  2048  azimuth  saii^les).  The  image  is  divided  into 
intervals  of  64  range  samples  and  2048  azimuth  sanqiles. 

2.  The  raw  SAR  data  area  is  range  compressed  and  comer  turned  yielding  a  dam  area  of  1024  R 
X  2048  A  sanqiles.  It  has  been  found  necessary  to  use  a  high  degree  of  inteipolation  (8  new 
samples  per  old  one)  to  reduce  the  modulation  effect  on  the  spectrum  caused  by  the  range 
sampling  grid.  Thus  the  input  to  the  azimuth  conqnession  stage  consists  of  (64  x  8)  R  x 
2048  A  samples. 
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3.  Based  on  the  ancillary  data,  initial  estimates  of  the  Doppler  rate  fpg  and  the  Doppler 
centroid  are  determined. 

4.  To  determine  the  range  dependence  of  the  range  compressed  data  are  divided  into  16 
range  intervals  (each  containing  64  x  8  range  lines,  where  8  is  due  to  the  interpolation).  An 
image  of  2048  azimuth  samples  is  obtained  using  the  full  available  ^leiture  of 4096  azimuth 
reference  samples.  For  each  interval  the  parameters  used  are  fpg  obtained  in  step  3  and  one 
of  16  values  of  fpp.  An  initial  value  of  fpp  is  obtained  in  step  3.  The  16  values  fpp  are 
chosen  to  be  equidistant  and  symmetrically  distributed  around  the  fpc  value  obtained  in  step 
3,  the  increment  is  S.SHz. 

5.  Azimuth  frequency  spectra  are  generated  by  obtaining  the  periodograms  of  the  complex 
imagery  in  the  azimuth  direction.  The  periodograms  are  averaged  to  give  a  higher  spectral 
fidelity.  An  example  of  spectra  for  different  processor  values  fpp  (point  2  in  Figure  1 .8)  is 
shown  in  Figure  1.3. 

6.  The  energy  difference 


A£  = 


is  computed  and  the  pair  A  E(Jpp),  fpp  is  stored.  An  example  of  a  A  £  versus  fpp  curve 
(point  3  in  Figure  1.8)  is  shown  in  Figure  1.9. 

7.  A  straight  line  is  fitted  to  the  16  pairs  (A  E,fpp),  and  the  fpp  value  corresponding  to  A  £  =  0 
(Jpc  i)  is  determined.  The  pair  (fp^  ,-4)  is  stored. 

8.  The  steps  4  through  7  are  repeated  for  all  16  range  intervals.  To  the  observed  (fpcpO  data,  a 
linear  function  is  fitted.  The  line  so  obtained  is  the  Doppler  line.  An  example  of  the  Doppler 
centroid  estimates  and  a  linear  fit  for  the  16  range  intervals  is  shown  in  Figure  1.10. 

To  get  an  estimate  of  the  error  in  each  determined  Doppler  centroid  value  the  approach  of  [Li  85] 
is  followed.  As  mentioned  before,  over  the  relatively  small  area  of  die  scene  of  ••  16km  by  8km, 
the  theoretical  variation  of  the  Doppler  centroid  as  a  function  of  range  is  approximately  linear 
(typical  nonlinear  effects  are  estimated  to  be  0.2Hz  in  the  reference).  So  the  root-mean-square 
deviation 


Raw  SAR  Data 


TNO  report 


Page 

25 


Figure  1.9:  Plot  of  energy  difference  against  several  values  of  processor  Doppler  centroid  fpp. 

Estimated  centroid  is 

from  the  regression  line  (the  Doppler  centroid  line  from  step  8)  is  taken  as  a  measure  of  accuracy. 
The  value  obtained  in  step  8  is  16.8Hz. 

A 

In  [Madsen  86],  the  standard  deviation  of  the  estimator /ipf  is  calculated  under  the  assumption 
that  the  scene  is  quasi  stationary,  that  is,  Ogfr*  ),  the  average  backscatter  coefficient  is  slowly 
varying  compared  with  the  system  inqtulse  reqtonse.  One  of  the  most  important  results  is  diat  the 
variance  of  the  estimator  is  proportional  to  the  scene  contrast,  SC,  defined  as 


SC  = 


<o“(r)>’ 


where  <  >  indicates  spatial  averaging. 
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Estiination  results  in  [Li  85]  and  [Curlander  82]  show  that  the  estimation  accuracy  is  highly 
dependent  on  scene  type,  e.g.  urban,  desert,  mountain,  etc.  and  that  successful  implementation 
requires  an  algorithm  capable  of  automated  detection  of  estimation  errors  that  result  from 
unsatisfactory  selection  of  a  test  area. 

Scenes  that  exhibit  an  asymmetric  azimuth  reflectivity  will  produce  a  degraded  Doppler  centroid 
estimate.  The  scene  types  that  exhibit  this  directional  reflectivity  characteristic  generally  have  a 
very  rough  surface  such  as  mountain  or  urban  areas.  Conversely  a  flat  surface  such  as  ocean  or 
desert  produce  highly  precise  estimates.  Test  results  in  [Curlander  82]  show  that  the  estimate 
uncertainty  for  ocean  is  between  l-3Hz  (1  a)  while  that  for  urban  areas  is  16-24//Z  (1  a). 


Doppler  centroid  (in  Hz) 


range 


Figure  1.10:  Doppler  centroid  estimates  and  linear  fit  for  16  range  intervals.  AE-tnethod. 
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4.3  Doppler  centroid  estiination,  time  domain  estimation 

An  alternative  way  of  finding  the  Doppler  centroid  is  by  using  the  Fourier  relation  between  the 
PSD  and  the  autocorrelation.  [Madsen  89]  describes  two  me±ods,  the  correlation  Doppler 
estimator  (CDE)  based  on  a  linear  estiniator  of  the  autocarrelation  function  and  the  sign-Doppler 
estiniator  (SDE)  based  on  a  nonlinear  estimator  of  the  autocorrelation  function.  The  SDE-method 
derives  the  autocorrelation  properties  of  an  assumed  circular  symmetric  Gaussian  process  by 
analyzing  the  sign  alone.  This  means  first  of  aU  that  the  majority  of  calculations  involved  in  the 
estimation  are  extremely  simple  (that  is,  sign  comparisons)  and  therefore  the  algorithm  is  well 
suited  to  real-time  implementation.  Secondly,  the  estimator  does  not  give  more  weight  to  strong 
targets  which  makes  it  less  sensitive  to  strongly  varying  scenes. 

The  conclusions  of  [Madsen  89]  are  that  the  time  domain  algorithms  are  efficiem  with  respect  to 
requirements  on  calculations  and  memory  and  hence  well  suited  to  real-time  systems  where  the 
Doppler  estimation  is  based  on  raw  data.  For  off-line  processors  where  the  Doppler  estimation  is 
performed  on  the  image  data,  the  A  £  estimator  and  the  CDE  algorithm  give  similar  performance. 
However,  when  it  comes  to  nonhomogeneous  scenes,  it  is  found  that  the  nonlinear  SDE  algorithm 
gives  superior  performance. 

We  have  implemented  both  the  CDE  and  the  SDE  algorithms.  Test  on  a  small  sample  of  data 
showed  results  almost  equal  to  the  ones  obtained  with  the  A  £-tnethod.  To  verify  the  claim  of  the 
SDE-algorithm  with  respect  to  nonhomogeneous  scenes  a  lot  of  additional  work  needs  to  be  done. 

4.4  Doppler  frequency  rate  determination,  AUTOFOCUS 

Sar  data,  in  the  azimuth  dimension,  have  the  property  that  from  any  subaperture  of  the  syrrthetic 
aperture  an  image  can  be  computed  of  the  same  target.  The  images  will  be  of  the  same  target 
except  that  the  resolution  will  be  reduced  because  only  a  part  of  the  full  ^rerture  (and  hence  from 
the  full  Doppler  induced  FM  bandwidth)  is  taken.  Use  of  an  incorrect  Doppler  frequency  rate 
value  fr,  in  the  processor  causes  the  defocus  of  the  image  and  a  diiff  in  the  azimuth  direction 
(also  called  map-drift).  This  shift  depends  on  the  position  of  the  subtqterture  in  the  synthetic 
aperture  (Section  1.3.2).  The  autofocus  method  is  based  on  this  shift  property.  The  aperture  is 
divided  into  two  subqiertutes  and  an  image  is  formed  from  each  subaperture  based  on  the 
processor  Doppler  frequency  rate 

The  two  inuiges  are  cross-correlated  to  detect  the  registration  error.  Based  on  the  detected  error  a 
new  processor  parameter  is  computed.  This  iteration  procedure  is  stopped  when  the  two  images 
become  focussed,  that  is  show  an  acceptable  diifi  with  respect  to  each  other.  This  mtq>-drifl 
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method  is  a  method  specific  to  measuring  a  quadratic  phase  error.  In  general,  n  sub^rertuies  are 
required  to  measure  an  nth-order  phase  error  polynomial  [Brown  88]. 

Figure  1.1 1  shows  the  block  diagram  of  the  Doppler  ftequency  rate  line  (the  Doppler  frequency 
rate  is  a  function  of  range)  estimator.  The  estimator  is  tested  with  Seasat  scene  data  from  an  area 
near  the  Goldstone  tracking  station.  The  scene  contains  several  bright  point-like  images  of  comer 
reflectors  that  were  deployed  for  calibration  purposes.  The  autofocus  algorithm  operates  as 
follows: 

1 .  Based  on  a  quick  look  processing,  a  portion  of  the  raw  S  AR  data  to  be  processed  into  an 

image  is  selected.  The  portion  is  divided  into  range  intervals,  the  data  in  an  interval  i  is  used 
to  estimate  the  average  cross-correlation  function  Over  the  relatively  small  area 

used  to  estimate  the  correlation  function,  the  Doppler  frequency  rats  /ojj  (which  is  a 
function  of  range,  formula  1.22,  is  considered  constant  fpg  (r)  s  fpjf  (i).  Since  the  processor 
rate  value  f^pit)  is  taken  constant  in  the  ith  range  interval,the  nuqj  drift  value  of  each  line  of 
the  interval  is  approximately  constant.  The  image  used  covers  an  area  «  37  km  by  8  km 
(2048  range  x  2048  azimuth  samples).  The  image  is  divided  into  128  intervals  of  16  lines  of 
2048  azimuth  samples. 

2.  The  raw  SAR  data  area  is  range  compressed,  with  1  x  interpolation  and  comer  turned 
yielding  a  data  area  of  2048  R  x  2048  A  complex  samples. 

3.  Based  on  the  ancillary  data,  initial  estitnates  of  the  Doppler  centroid  fpc  l^e  Doppler 
frequency  rate  fpp  are  determined. 

4.  A  range  interval  i  is  azimuth  compressed  into  an  L-look  and  an  R-look  image  of  16  x  2048 
samples.  The  complete  available  spectral  band  of  4096  azimuth  samples  is  used.  The  band  is 
divided  into  two  non-overlapping  subbands  of  2048  azimuth  sanqtles.  The  processor 
parameters  used  are  fpc  as  obtained  in  step  3  and  the  interval  Doppler  frequency  rate  value 
frpii)  =  A  +  Bi.  An  initial  value  of  f^p  is,  ffp{i)  s  f^p  (A=fpp  as  obtained  in  step  3  and 
B=0).  Figure  1.12  shows  an  example  of  look  images  obtained  by  using  three  different 
processor  values  for  the  Doppler  frequency  rate.  Here  the  range  interval  consists  of  80 
samples.  For  each  line  the  same  fpp  value  (Doppler  centroid  value  used  by  the  processor)  is 
taken. 

5.  The  L-look  and  R-look  images  are  line  by  line  cross-correlated  and  the  obtained  16  cross¬ 
correlation  functions  are  averaged  yielding  Ru^kJ).  Typical  average  cross-correlation 
functions  R^(k)  for  different  values  of  the  Doppler  frequency  rate  used  by  the  processor  are 
shown  by  Figure  1.13.  Due  to  the  blurring  of  the  look-images  when  incorrect  fpp  values  are 
used,  the  cross-correlation  peak  broadens. 
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6.  The  location  kmax(C)  of  the  peak  of  the  average  cross-conelation  function  Ru({k,i)  (i.e.  the 
shift  between  both  images)  is  determined.  The  quality  of  the  lag  estimate  deteriorates  as  the 
error  in  the  processor  Doppler  frequency  rate  increases  due  to  the  broadening  of  the  cross- 
correlation  peak. 


I  . 
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0  l-look 
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0  R-look 
Image 


in  Average 
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^  Observed  Lag  Line 


^  Estimated  Lag  Line 


Doppler  Rate  Line 
fDR{i)  =  A  +  6i 


Figure  1.11; 


Block  diagram  of  Doppler  rate  estimation  algorithm  (AUTOPOCUS) 
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7.  The  steps  4  through  6  are  tepeaxed  till  all  128  range  intervals  have  been  processed.  To  the 
observed  data  kiiiax(i),  a  linear  function  a  -«■  b  i  is  fitted.  In  the  ideal  case,  perfect  focussing 
occurs  when  kmax  >0,  a  =  0.b  =  0.  Inthe  beginning  of  the  iteration  procedure,  in  general,  a 
and  b  will  substantially  deviate  from  zero.  The  determined  values  then  are  used  to  compute 
new  processor  constants  A  and  B  to  be  used  in  =  A  -f  B  i,  and  a  next  iteration  is 
started.  Figure  1.14  shows  lag  lines  observed  for  different  processor  Doppler  fiequency 
functions,  the  top  and  bottom  functions  correspond  with  heavily  defocnssed  images,  the 
middle  one  ccrresponds  with  the  focussed  images.  The  iteration  procedure  comes  to  an  end 
wheii  the  stop  criterion  £^kmax(i)2  <  e  is  met  The  corresponding  processor  values,  A 
and  B,  are  the  best  estimates  of  the  Doppler  rate  line/^j^fi) »  A  -t-  B  i.  The  lag  estimation 

error  decreases  when  the  width  of  the  cross-correlation  function  decreases.  A  sharp  peak  will 
give  smaller  errors  than  a  wide  one.  Thus,  the  error  is  not  only  related  to  the  signal-to-noise 
ratio,  but  also  to  the  scene  characteristics.  The  autofocus  technique  is  most  accurate  for  a 
scene  that  exhibits  bright  point  targets  and  is  significantly  degraded  for  ocean  or  desert.  The 
scenes  with  relatively  few  bright  targets  will  show  broader  cross-cmrelation  peaks  and 
therefore  have  a  higher  degree  of  uncertainty  for  the  Doppler  fiequency  rale  estimate. 

In  [Li  8S]  test  results  are  shown  for  the  autofocus  technique  for  the  Goldstone  tracking 
station  scene,  an  agricultural  field  pattern  and  a  patch  of  ocean  in  the  Bay  area.  The  rms 
deviations  are  «•  0.1  Hz/sec  for  the  Goldstone  scene  and  0.2  Hz/sec  for  the  agricultural  scene. 
The  increased  deviafitm  in  case  of  the  agricultural  scene  is  attributed  to  the  lack  of  bright 
point  Uke  targets  in  the  scene.  Visual  examination  of  the  ocean  image  revealed  no  sharp 
boundaries  or  point-like  features,  and,  no  consistent  trend  in  the  observed  maxima  was 
found.  They  concluded  that  the  autofocus  method  cannot  be  applied  to  this  type  of  scene.  For 
most  images  ovct  land,  however,  they  believed  that  this  method  Aould  be  effective. 

The  rms  deviation  of  our  processing  result  was  0.04  Hz/sec. 


kmax+51  kmax+51 


range  interval 


Fiyure  1.14: 


range  interval 


Obaerved  laa  lines  for  diffetcMnoccMorDnipla' frequency  liBic«ioat/^i)B  A -f  B  L 
From  top  to  tioaom  Ai^SiK ,  B«D ;  A>S18 .  B^.03  mlamty  fbcotied;  A^26 ,  B^. 
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A.  1  RANGE- AZIMUTH  SHIFT  INDUCED  BY  AN  INCORRECT  DOPPLER 

CENTROID  VALUE  OF  THE  PROCESSOR. 

Figure  A.l  shows  the  position  of  the  point  target  locus,  t  =  0  corresponds  with  dre  moment  where 
the  target  is  in  the  antenna  boresight  direction.  Also  shown  is  the  reference  locus  with  centre 
(-t.p). 

At  the  time  t  an  image  point  (in  general  not  the  image  point  of  the  point  target)  is  computed  by 
the  azimuth  compression. 

The  correlator  output  due  to  the  point  target  is 


C  ( T )  =  s,  ( T -I- p  )s  J  ( p  )dp 


(A.l) 


where  is  the  azimuth  integration  time.  s,(P)  the  azimuth  reference  function  and  +  P)  the 
point  target  signal  encountered  along  the  reference  locus.  If  the  processor  uses  the  correct  Doppler 
parameters,  the  Doppler  centroid  value  f/y-  and  Doppler  rate  the  maximum  point  target 
output  occurs  at  T  =  0  (and  p  =  0).  If  incorrect  processor  parameter  values  are  used,  the  maximum 
point  target  response  will  occur  at  a  value  of  t  ^  0  (and  p  *■  0).  The  value  of  T  and  p  as  a  function 
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of  parameter  mismatch  is  derived  below. 

The  target  locus  is 

r,(t)  =  ao, +a„t  +  a2,t^ 

The  reference  locus  is 

r,(P)  =  r,(T)  +  aifP  +  a2rP^ 

For  the  difference  between  target  and  reference  locus  we  find 

8,  (P)  =  r,  (T  +  P)  -  r^(P)  =  r,(T  +  P)  -  r,(T)  - 1  fl,,P  +  a2rf ) 

From  A.2  it  follows 

r,  (t  +  P)  =  r,  (T)  +  (ai,  +  2a2,T)p  +  32, P^ 

Subsumtion  of  A.5  into  A.4  gives 

5,  ( P )  =  (a,,  +  2  a2,T  )P  +  a2,p^  -  (aj^P  +  a2,P^ ) 

Assuming  a  correct  processor  value  for  <22  we  have 

32,  =  32, 

hence 

8r(P)=((ai,  +2a2,T)-ai,)P 

We  see  that  the  loci  coincide  for 

-  _£lLZilL 
232, 

With 

j,(a)  =  a(a)exp7<t»,(a) 

where  <i(a)  is  caused  by  the  antenna  and 

j,(P)  =  exp  7(1),  (P) 

with 

(t>,(a)  =  Co,  +c„a  +  Cj,a* 

0r(P)  =  f<>  +c„P  +  Cj,p’ 


(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 
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we  find 

C(T)  =  exp  j(co,  +  Ci,T  +  C2tX^  -  Cor ) 

from  1.7,  .29  and  .31 


^2,  = 

it  follows 


C(to)  =  e!tpj{(ti,(Xo)-CorljJ*^^^a(to+P)<iP  (A. 11) 

Thus  the  coincidence  of  the  target-  and  reference  locus  yields  an  "inphase"  integratioo  and  thus 
the  maximum  output. 

Note  1 :  We  tacitly  assume  the  “inphase"  integration  to  be  far  more  inqxntant  than  a  value  of  T 
smaller  than  x,.  It  is  true  that 

<»  ( 7  +  P )  d  p  <  a(x<,  +  P)  d  P 

when 

x<x„ 

but  since  the  reference  locus  and  the  target  locus  do  not  coincide,  a  jdiase  function  occurs  and  the 
signal  amplitude  along  the  reference  locus  is  not  coherently  added. 

IC(T)l=l£‘%xp7[e(T)P]a(T  +  P)dpi<£'”a(T,-t-P)dP)l 

Note  2:  From  the  form  of  the  antenna  pattern  and  .33  follows 

IC(Xo)l<  a(P)dp  (A. 12) 

i«/2 

so  due  to  the  incorrect  a,,  value  used  we  lose  signal  energy. 

The  range  offset  value  p  which  ctmesptxtds  to  x,  is 


TNO  report 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE  (MOD-NL) 


1 .  DEFENSE  REPORT  NUMBER  (MOO-NL)  2.  RECIPIENTS  ACCESSION  NUMBER 

TD91-3978 


3.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

FEL-91-B331 


4.  PROJECT/TASKANORK  UNIT  NO.  5.  CONTRACT  NUMBER 

22020 


6.  REPORT  DATE 
NOVEMBER  1991 


7.  NUMBER  OF  PAGES 

39  (INCL  APPENDIX,  EXCL.  RDP  & 
DISTRIBUTION  LIST) 


8.  NUMBER  OF  REFERENCES 
8 


9.  TYPE  OF  REPORT  AND  DATES  COVERED 
FINAL 


10.  TITLE  AND  SUBTITLE 

ESTIMATION  OF  THE  PARAMETERS  USED  IN  THE  AZIMUTH  COMPRESSION 


11.  AUTHOR(S) 

F.P  PH.  DE  VRIES 


12.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

TNO  PHYSICS  AND  ELECTRONICS  LABORATORY,  P.O.  BOX  96864,  2509  JG  THE  HAGUE 
OUDE  WAALSDORPERWEG  63.  THE  HAGUE,  THE  NETHERLANDS 


13.  SPONSORINQimONITORING  agency  NAME(S) 

NETHERLANDS  AGENCY  FOR  AEROSPACE  PROGRAMS  (NIVR) 


14.  SUPPLEMENTARY  NOTES 


15.  ABSTRACT  (MAXIMUM  200  WORDS,  1044  POSITIONS) 

IN  AZIMUTH  PROCESSING  OF  SAR  DATA,  THE  SLANT  RANGE  HISTORIES  OF  THE  TARGETS  ARE  REQUIRED.  AFTER  RANGE 
COMPRESSION  THE  MAXIMA  OF  THE  RESPONSES  OF  A  POINT  TARGET  ARE  ON  THE  SO  CALLED  RANGE  MIGRATION 
LOCUS  FOR  SPACEBORNE  SAR  THE  LOCUS  TRAVERSES  A  NUMBER  OF  RANGE  BINS.  THIS  LOCUS  DETERMINES  WHERE 
TO  GET  THE  SAMPLES  NEEDED  FOR  THE  AZIMUTH  COMPRESSION  AND  THE  PHASE  FUNCTION  WHICH  MUST  BE  USED. 
FOR  SATELLITE  SAR  THE  LOCUS  CAN  BE  DESCRIBED  BY  TWO  PARAMETERS.  THE  DOPPLER  CENTROID  fpc  AND  THE 
FREQUENCY  RATE  THESE  PARAMETERS  CAN  BE  ESTIMATED  BY  USING  ACCURATE  SPACECRAR  ANCILLARY  DATA 
AND  BY  UTILIZING  THE  RECEIVED  RADAR  DATA  ALONE.  THIS  REPORT  DESCRIBES  THE  EFFECTS  OF  ERRORS  IN  THE 
PARAMETERS  USED  BY  THE  PROCESSOR  AND  DISCUSSES  ALGORITHMS  TO  ESTIMATE  THE  PARAMETERS.  EXAMPLES  OF 
RESULTS  FROM  PROCESSING  SEASAT  DATA  ARE  SHOWN. 

THE  REPORT  IS  ONE  OF  A  SERIES  FOR  A  SAR  PROCESSOR  FOR  THE  ERS-1  DATA,  PHASE  2,  PRECISION  PROCESSOR.  THE 
WORK  DESCRIBED  WAS  SUBSIDIZED  BY  THE  NIVR  (NIVR.  NRT  NR,  2802  FE.  SAR  PRIXIESSOR). 


16.  DESCRIPTORS 

SYNTHETIC  APERTURE  RADAR 
SATELLITE  BORNE 
ESTIMATION 


IDENTIFIERS 

DIGITAL  SIGNAL  PROCESSOR 
DOPPLER  PARAMETERS 
DOPPLER  CENTROID 
DOPPLER  FREQUENCY  RATE 


17b.  security  CLASSIFICATION 
(OF  REPORT) 

UNCLASSIFIED 


17b.  security  CLASSIFICATION 
(OF  PAGE) 

UNCLASSIFIED 


17c.  SECURITY  CLASSIFICATION 
(OF  ABSTRACT) 

UNCLASSIFIED 


18.  DISTRIBUnON/ AVAILABILITY  STATEMENT 
UNLIMITED 


17(1  SECURITY  CLASSIFICATION 
(OF  TITLES) 

UNCLASSIHED 


